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Abstract 

We present a numerical analysis of genus statistics for dark matter halo catalogs from the Hubble 
volume simulation. The huge box-size of the Hubble volume simulation enables us to carry out a reliable 
statistical analysis of the biasing properties of halos at a Gaussian smoothing scale of Rq > 30/i~^Mpc 
with a cluster-mass scale of between 7 x 10^^h~^ Mq and 6 x lO^'^/i~^M0. A detailed comparison of the 
genus for dark matter halos with that for the mass distribution shows that the non-Gaussianity induced 
by the halo biasing is comparable to that by nonlinear gravitational evolution, and both the shape and the 
amplitude of the genus are almost insensitive to the halo mass at Rq > 30/i~^Mpc. In order to characterize 
the biasing effect on the genus, we apply a perturbative formula developed by Matsubara (1994). We 
find that the perturbative formula well describes the simulated halo genus at Rq > 50/i~^Mpc. The result 
indicates that the biasing effect on the halo genus is well approximated by nonlinear deterministic biasing 
up to the second-order term in the mass density fluctuation. The two parameters describing the linear and 
quadratic terms in the nonlinear biasing accurately specify the genus for galaxy clusters. 

Key words: cosmology: dark matter — cosmology: large-scale structure of universe — galaxies: 
clusters: general — galaxies: halos — methods: statistical 



1. Introduction 

The genus statistics characterizes the topological na- 
ture of a density field. In contrast to the more conven- 
tional two-point statistics, such as the two-point correla- 
tion functions or the power spectrum, the genus is sen- 
sitive to the phase information of a density field; there- 
fore, the genus provides complementary information of 
the present cosmic structure. The genus has an analytic 
formula in a density field described by random-Gaussian 
statistics. Since Gott et al. (1986) proposed the genus 
of the large-scale structure as a test of the hypothesis 
that the initial perturbations are Gaussian, the genus has 
been evaluated for various galaxy or galaxy cluster cat- 
alogs (e.g., Gott et al. 1989; Park et al. 1992; Moore 
et al. 1992; Rhoads et al. 1994; Vogeley et al. 1994; 
Canavezes et al. 1998; Hoyle et al. 2002; Hikage et al. 
2002). These investigations generally indicate consistency 
with the initial Gaussianity. 

For a statistical analysis based on the distribution 
of luminous objects, however, non-Gaussianity produced 
in the process of cosmic structure formation should be 
properly considered. Gravitational nonlinear evolution 
is one of the most important non-Gaussian sources. It 
has been well studied using a variety of models, in- 
cluding a perturbative correction by the Edgeworth ex- 
pansion (Matsubara 1994), direct numerical simulations 
(Matsubara, Suto 1996), an empirical lognormal model 
(Matsubara, Yokoyama 1996) and the Zcl'dovich approx- 
imation (Seto et al. 1997). 

Another non-Gaussian source is biasing, i.e., the sta- 



tistical uncertainty between the spatial distribution of lu- 
minous objects and that of dark matter. Biasing is gen- 
erally described by a nonlinear and stochastic function 
of the underlying mass density from various recent work 
(Dekel, Lahav 1999; Blanton et al. 1999; Somervillc ct 
al. 2001; Taruya et al. 2001). The difficulty in a gen- 
eral prediction of a biasing expression comes from the 
fact that the biasing properties are sensitive to the com- 
plicated formation mechanism of luminous objects. On 
the other hand, the biasing for dark matter halos, which 
are the most likely sites for galaxies and clusters, is rela- 
tively easy to be modeled and thus well investigated on the 
basis of both analytical and numerical approaches (Mo, 
White 1996; Taruya, Suto 2000; Yoshikawa et al. 2001). 
In particular, dark halos with cluster-scale mass are usu- 
ally supposed to have a fairly good one-to-one corre- 
spondence with galaxy clusters, at least in a statistical 
sense. Therefore, biasing models for dark halos can be 
applied to upcoming galaxy cluster catalogs with an em- 
pirical relation between the halo mass and the tempera- 
ture/luminosity of galaxy clusters (e.g., Kitayama, Suto 
1997). In this spirit, Hikage et al. (2001) derives an an- 
alytical model of genus statistics for galaxy clusters with 
a nonlinear and stochastic halo biasing model of Taruya 
and Suto (2000). 

A complementary approach to halo biasing is a direct 
analysis of cosmological numerical simulations. The sta- 
tistical significance of previous studies, however, has been 
limited by the number of particles and the box-size cm- 
ployed in the simulations. Since the mean separation of 
rich clusters of galaxies is ~ 50/i^^Mpc, a typical cosmo- 
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logical simulation in a SOO/i^^Mpc box has merely - 200 
clusters. With this number of clusters, it is almost im- 
possible to detect the cluster-mass dependence of various 
statistical measures in a reliable manner (e.g., Watanabe 
et al. 1994). This situation may be significantly im- 
proved by using the Hubble Volume Simulation (Jenkins 
et al. 2001). Specifically, we use the ACDM (Lambda Cold 
Dark Matter) data (AHVS, hereafter) at z = 0, which em- 
ploy N = 10^ particles in a box-size of 3/i~^Gpc. The data 
enable us to quantitatively address the mass and scale de- 
pendence of the halo biasing and the related statistics. 

In the present work, we perform a numerical analysis of 
the genus for dark halos with the Hubble volume simula- 
tion in order to clarify the uncertainty of the halo biasing 
in detecting the initial non-Gaussianity. From a detailed 
comparison between the genus for dark halos and that 
for dark matter, we find that the halo biasing effect on 
the non-Gaussianity of genus is comparable to the non- 
Gaussianity induced by the nonlinear gravitational evo- 
lution. In order to characterize the halo biasing effect 
on the genus statistics, we apply a perturbative formula 
developed by Matsubara (1994), and compare it with 
the simulated genus. We find that the perturbative for- 
mula properly describes the halo genus at a smoothing 
scale larger than 50ft.~"'^Mpc. Also, the biasing effect on 
the genus for halos is shown to be well-approximated by 
the nonlinear deterministic biasing up to the second-order 
terms in the mass density fluctuation. 

The rest of the paper is organized as follows. We briefly 
describe the simulation data that we use and the halo 
catalogue in section 2. We then discuss in section 3 the 
biasing properties of dark halos while paying particular 
attention to their nonlinear and stochastic nature. We 
present a detailed study of the topology of the spatial 
distribution of dark halos using the genus statistics based 
on numerical simulations and a perturbative analysis in 
section 4. Finally section 5 is devoted to conclusions. 

2. Dark Matter halos in the Hubble Volume 
Simulation 

In the present work, we use the z — snapshot data 
of ACDM. The model assumes the density parameter 
of r^o = 0.3, cosmological constant of Xq = 0.7, and the 
Hubble constant in units of 100km s~^Mpc~^, h = 0.7. 
The fluctuation amplitude, smoothed over a top-hat ra- 
dius of 8/i~^Mpc, is set as erg = 0.9 so that the model 
satisfies both the cluster abundance and the COBE nor- 
malization (e.g., Eke et al. 1996, Kitayama, Suto 1997). 
The initial density field is computed with CMBFAST 
assuming baryon density parameter of ^l\^h^ = 0.0196, 
and the random-Gaussian statistics (Colberg et al. 2000). 
The simulation box-size is 3000/i~^Mpc and the softening 
length for the gravitational force is set to be 100/i^^kpc. 
An initial distribution of = 10^ dark matter particles 
is generated using the Zel'dovich approximation, and is 
advanced with the Particle-Particle-Particle-Mesh (P'^M) 
code. The mass of each dark matter particle is set at 
2.25 X 10^^h~^MQ. Hamana et al. (2001) made sure that 



the two-point clustering feature at z = in the AHVS is 
reliable even at an order of magnitude below the mean 
particle separation 3h~^Mpc (still a few times larger than 
the adopted softening length). Thus, the dynamic range 
of the simulation data amounts to ~ 10^ in length. 

Dark halos are identified with a friends-of-friends finder 
with a linking length of = 0.164 in units of the mean- 
particle separation. We define linked groups with the min- 
imum number of member particles to be 30 as dark ha- 
los in the current analysis. We divide those halos (about 
1.5 X 10^ in total) in five different subgroups according to 
their mass (table 1) so that each subgroup includes about 
3 X 10^ halos. Because of the huge size of the AHVS, 
the number of halos in each subgroup is sufficiently large 
to characterize the mass dependence of various statistical 
measures in a reliable manner. In the analysis discussed 
below, we evaluate the error for each statistical measure 
from the sample-to-sample variation among the results for 
eight subboxes (box-size of 1500ft.~^Mpc) constructed out 
of the entire simulation box. 



Table 1. List of five subgroups of dark matter halos in the 
AHVS. 



subgroup 


mass ranj 




number of halos 


SS 


6.7 X 10'^ 


- 8.0 X lO^'^ 


328449 


S 


8.0 X 10^3 


- 1.0 X 10" 


321732 


M 


1.0 X 10^^ 


- 1.3 X 10" 


295204 


L 


1.3 X 10" 


- 2.0 X 10" 


300828 


LL 


2.0 X 10" 


- 6.0 X lO^^ 


314782 



3. The Nature of Biasing of Dark Halos 

First we consider the nonlinear and stochastic nature 
of halo biasing with particular attention paid to its mass- 
and scale-dependence. For this purpose, we compute the 
joint probability distribution function (PDF) of the num- 
ber density field of halos, (5haio, and the mass density field 
of dark matter particles, (^mass- Those density fields are 
first evaluated at 128'^ grids for each subbox data of par- 
ticles smoothed over a Gaussian radius of Rq. Then, the 
averaged value of the joint PDF over eight subboxes is 
computed. Figure 1 shows the averaged joint PDF for 
specific halo subgroups, LL and SS, at smoothing scales 
of 30 and bOh'^Mpc. 

The dotted lines represent the mean biasing, 
<5haio('5mass), averaged over 5haio at a fixed value of 

^mass • 

<^halo('5mass) = J (^halo ^'('^halo If^mass) d^halo- (1) 

Here, /'((5haio|i5mass) represents the conditional PDF of the 
halo fields at a given ^mass, which is equal to the joint 
PDF, P(5mass, <5haio), divided by the one-point PDF of 
dark matter, P((5mass)- We fit the mean biasing to the 
following quadratic nonlinear model: 

^bias,fit(^mass) — ^If^niass H (<^mass f^mm ); (2) 
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Fig. 1. Joint probability distribution functions, 
f'('5mass, 5halo)i o{ the halo and the dark matter den- 
sity fields computed from the AHVS data with Gaussian 
smoothing; Upper-left: LL halos with Rq = 30/i~^Mpc, 
Upper-right: LL halos with Rq = 50/i~^Mpc, Lower-left: 
SS halos with Rq = 30h~^Mpc, Lower-right: SS halos with 
Rq = 50/i~^Mpc. Contour lines show P {Smaee , Si^alo) = 10~^, 
10~^ and 10~^. Dotted, dashed and solid lines indicate the 
mean biasing [equation (1)] and the linear and the quadratic 
fits [equation (2)] to the mean biasing respectively. 

where crmm = (<^mass^)^^^ and (• • •) denotes the average over 
P{Smass)- The fitted values for the hnear and the second- 
order biasing coefficients, bi and 62, in each subgroup are 
Hsted in table 2. While a linear fit (dashed lines) with 
62 = in equation (2) fails to approximate the simula- 
tion result, a quadratic fit (solid lines) is quite acceptable 
for our interest range in the biasing of cluster-scale halos 
(Hikage et al. 200f). 

On the other hand, the stochasticity in the biasing is 
significantly affected by a shot-noise term in these halo 
samples due to the small number density of halos. For ex- 
ample, the number of halos within the Gaussian smooth- 
ing volume, (27r)'^/^i?Q is around 22 at Rq ~ 50h~^Mpc, 
and then the shot-noise term amounts to 0.21 (= 1/^/22), 
which is comparable to the r.m.s. variance listed in ta- 
ble 2. ft is, therefore, not easy to extract the information 
concerning the pure stochasticity in the clustering biasing 
from our halo catalogs. We therefore focus on the non- 
linear biasing effect and compare the genus for dark halos 
with perturbativc predictions. The effect of the stochas- 
ticity of biasing on the halo genus is discussed in appendix 
1 using a perturbativc analysis. 

4. Genus Statistics for Dark Halos in AHVS 

4-1. Definition and Computing Method of Genus 

Genus, g{vcr), is defined as —1/2 times the Eulcr char- 
acteristic of the isodensity contour of the density field 5 



at the threshold level of Va times the r.m.s. fluctuation 
(T. In practice, this is equal to (number of holes) — (num- 
ber of isolated regions) of the isodensity surface. In the 
random-Gaussian field, genus per unit volume is given by 
the following analytical form: 

where cri ee (iVJp)^/^ and a = (J^)^/^ {V5 denotes the 
spatial derivative of the density field 5 and (• • •) denotes 
the average over of the PDF of 5 and V5 (cf., Doroshkevich 
1970; Adler 1981; Bardeen et al. 1986; Hamilton et al. 
1986). 

In contrast, the genus is often expressed as a function of 
Vi instead of (e.g., Gott et al. 1989). The threshold ft 
is used to parameterize the fraction / of the volume that 
lies on the high-density side of the contour, 

f = -j= / e-^'/^dx. (4) 
V 27r j^f 

Figure 2 shows a comparison between and Vf. The devi- 
ation from a linear relation represents the non-Gaussianity 
of the one-point PDF in the halo field smoothed with 
scales of 30, 50 and 100/i"^Mpc from the left to right pan- 
els. While the genus with Rq = 100ft-~^Mpc approaches 
the random-Gaussian prediction, the non-Gaussianity of 
the one-point PDF is detectable with a smaller Rq and a 
larger mass of halo. Up to the first order of perturbativc 
correction in ct, the relation between and i^f is given by 
(Matsubara 2003) 

Va=Vf+—^{vi - 1)(T, (5) 

where we directly compute a skewness parameter, S''-''-' (= 
{v„^)/(j), from the halo field. Figure 2 shows that the 
lowest order approximation describes the simulated results 
very accurately at scales of our interest. 

If the evolved density field has an exact one-to-one cor- 
respondence with the initial random-Gaussian field, then 
this transformation from v„ to removes the effect of evo- 
lution of the one-point PDF of the density field. Under 
this assumption, the genus as a function of the volume 
fraction, expressed as <?(ff), is sensitive only to the topol- 
ogy of the isodensity contours, rather than evolution with 
time of the density threshold assigned to a contour. The 
limitations of the one-to-one mapping between the initial 
and evolved density fields are critically examined by Kayo 
et al. (2001). In the following section, wc present the 
results of the genus in terms of both and v^. 

The computation of genus proceeds as follows: we first 
assign densities of halos and dark matter particles on 128^ 
grids using the cloud-in-cell method. The distance be- 
tween each grid is sufficiently small so as not to affect the 
density field smoothed over a scale of 30/i~^Mpc. The 
density fields are Fourier-transformed, multiplied by the 
Gaussian window with a smoothing length of Rq, and 
then transformed back to a real space. The smoothed den- 
sity fields are used in defining the isodensity surface with a 
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Table 2. Linear and second-order biasing coefficients. 

















100/i"%Ipc 




subgroup 






K 


4 




K 


4 





"SS Tm -0.72(±0.05) (l32 L66 -0.69(±0.13) OlB TGS 1.50(±1.18) 006 

S 1.76 -0.60(±0.04) 0.33 1.76 -0.76(±0.14) 0.16 1.72 -0.24(±1.18) 0.06 

M 1.88 -0.31(±0.04) 0.35 1.87 -0.53(±0.15) 0.17 1.88 -1.95(±1.05) 0.06 

L 2.16 0.38(±0.05) 0.36 2.12 -0.27(±0.16) 0.18 2.05 -0.40(±1.08) 0.06 

LL 2.94 4.37(±0.03) 0.42 2.97 4.48(±0.11) 0.21 3.05 1.59(±0.96) 0.07 



* Linear biasing coefficients. 
^ Second-order biasing coefficients witli la fitting error [equation (2)]. 
^ the amplitude of tlie rms density fluctuation of lialos a = ((5haio^)^^^ with Gaussian smoothing. 



3 
2 
1 







I I I I I I I I I I I I I I I I I I 
i?c=30h-'Mpc 

LL.sim. 
; LL,eq.[5] 




I I I I I I I I I I I I I I I I I I I \j^y 

/?G=50/i-iMpc 



LL.sim. 
LL,eq.[5] 
SS.sim. 
SS,eq.[5]. 



I I I I I I I I I I I I I I I I I I I I I . 

LL.sim. 

: LL,eq.[5] 

SS.sim. 
SS,eq.[ 



-1 

-2 



~ a SS.sim. 
SS.eq. 



I I I I I I I I I I I I I I I I I I I I I I 




l'„=l': 




l'„ = l^, 



■1 







1 











Fig. 2. Comparison between U{ and directly computed from the density fields of LL (filled circles) and SS (open triangles) halo 
subgroups with sample-to-sample variations. The perturbative expression up to the first order of a [equation (5)] is also plotted 
for the LL subgroup (solid lines) and the SS subgroup(dotted lines). From left to right, Gaussian smoothing scale is 30,50, and 
100/i~^Mpc, respectively. The dashed lines indicate the linear relation. 



given threshold, and the resulting genus is evaluated with 
the CONTOUR 3D code (Weinberg 1988 ) by integrating 
the deficit angles over the isodensity surface. 

4-. 2. Genus Statistics for Dark Matter and Dark Halos 

Since the one-point PDF of dark matter particles is em- 
pirically known to be well approximated by the log-normal 
PDF (e.g. Coles, Jones 1991; Kofman et al. 1994; Kayo et 
al. 2001), it is natural to expect that it is also the case for 
genus. In fact, Matsubara and Yokoyama (1996) derived 
a genus expression for dark matter assuming the nonlinear 
density field of dark matter has a one-to-one mapping to 
its primordial Gaussian field. If one adopts the log-normal 
mapping, their result can be explicitly written as 



where the maximum value, 



5ln(!^(t) = .gMAX.LN[l - Xl^{l'„)] eXp 



,K) 



^ ln[(l^H^^aVl+^ 



(6) 



(7) 



1 



,^MAX,LN — 



(2^) 



is defined by 

3/2 



(8) 



_3(l + a2)ln(l + cr2) 

Note that in terms of just reduces to the Gaussian 
prediction (3). In reality, this one-to-one mapping as- 
sumption is much stronger than the statement that the 
one-point PDF is empirically well approximated by the 
log-normal PDF. 

We compare the genus for dark matter averaged over 
eight sub-boxes with the analytical expressions [equations 
(3) and (6)] multiplied by the volume of the sub-boxes, 
(1500/i~^Mpc)3. The calculations of cr and CTi in the an- 
alytical predictions are based on the power spectrum by 
Peacock and Dodds (1996). Figure 3 shows that genus as 
a function of [panel (a)], the difference, between the 
genus for the simulated dark matter ^dm and the analyti- 
cal predictions (/theory, normalized by the maximum value 
of ffDM, 5max,dm [panel (b)], and the same as panels (a) 
and (b), but as a function of [panels (c) and (d)]. 
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Fig. 3. Panel (a): Genus for dark matter computed from AHVS with filled circles and the analytical genus expression in ran- 
dom-Gaussian with dotted lines [equation (3)] and in log-normal statistics with solid lines [equation (6)]. The top, middle, and 
bottom panels, respectively, indicate the results for smoothing lengths of 30/i~^Mpc, 50/i~^Mpc, and 100h~^Mpc. This order of 
the smoothing scale is the same in the other panels. Panel (b): The difference between our simulation results (denoted by Qdm) and 
the theoretical predictions (denoted by gthoory; the random-Gaussian expression with dotted lines, the log-normal model with solid 
lines containing error-bars) normalized by the amplitude of the simulated genus (denoted by Smax.dm)- Panels (c) and (d): Same as 
panels (a) and (b), but for genus as a function of i'{ [equation (4)]. The simulated results arc compared with the random-Gaussian 
expression with error-bars (denoted by gRc) 




Fig. 4. Same as figure 3 but for dark halos in LL (filled circles) and SS (open triangles) subgroups. The sam- 
ple-to-sample variance of genus is plotted with error-bars only for the LL subgroup. Panels (a) and (b): Comparison 
between the genus for halos (denoted by gHaio) with the simulated genus for dark matter (denoted by ^dm) as a func- 
tion of UcT- Panels (c) and (d): Comparison between genus for halos as a function of U{ with equation (3). 
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We estimate the sample-to-sample variance of the sim- 
ulated genus Ag(i/) by 

(-8 ^1/2 

A5H = |E[-9'M--9H]V7| , (9) 

where gi{i') denotes the value of the genus on the contour 
surface at the threshold estimated from sub-box i; g{i') 
is the value of the genus averaged over eight sub-boxes. 

Panels (a) and (b) show that the empirical prediction 
(6) well explains our current results for dark matter in- 
cluding both the amplitude and the shape of the genus in 
our interested scale. This is consistent with the numeri- 
cal results by Matsubara and Yokoyama (1996). Panels (c) 
and (d) show that the genus for dark matter as a function 
of Uf agrees with the analytical prediction for a Gaussian 
random field [equation (3)] . 

We then plot the genus for dark halos of two specific 
subgroups, LL and SS in figure 4. While panels (a) and 
(b) show comparisons of the halo genus with the simu- 
lated genus for mass as a function of I'^r, the panels (c) 
and (d) represent the results in terms of i/f. The differ- 
ences in the genus for halos from the simulated genus for 
dark matter [panel (b)] or the random-Gaussian predic- 
tion [panel (d)] are normalized by the maximum value of 
the genus for halos, ^MAx.Haio. Figure 4 shows that the 
non-Gaussianity induced by halo biasing is comparable to 
that by the non-linear gravitational evolution, and that 
the shape and the amplitude of the halo genus are almost 
insensitive to the halo mass. The result suggests that the 
infiuence of biasing can be corrected perturbatively. To 
see the weak dependence of the genus on the halo mass 
more clearly, we plot in figure 5 the ratio of the genus am- 
plitude, gMAx,Haio/<?MAx,DM, i-G., tlic maxlmum value of the 
genus amplitude for dark halos divided by that for dark 
matter. The results for each halo subgroup arc plotted as 
a function of the averaged halo mass. While the increase of 
the smoothing scale significantly changes the genus ampli- 
tude in figure 4, the relative amplitude in figure 5 almost 
remains unchanged for a wide range of the halo mass. 
Furthermore, the amplitude of halo genus is very close to 
that of dark matter, with only 6-16 percent deviations 
at Rq ~ 30/i~^Mpc. This is indeed consistent with the 
prediction of Hikage et al. (2001) based on the one-to- 
one mapping assumption from the genus for dark matter. 
Therefore, one can conclude that the biasing effect on halo 
genus on the cluster-mass scale is generally small, and that 
the perturbative analysis based on Edgeworth formula is 
safely applicable at Rq > 30/i~^Mpc. 

4-3. Comparison with Perturbation Theory 

Having established the weak infiuence of the biasing, 
we wish to further examine the halo genus for a better 
understanding of the nature of biasing. For this purpose, 
we attempt to model the halo genus using a perturba- 
tive approach with particular emphasis on the nonlinearity 
of the biasing effect. We start with Matsubara's expres- 
sion for the genus density in a weakly non-Gaussian field 
(Matsubara 1994, 2003). In the lowest-order approxima- 
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Fig. 5. Ratio of the amplitude of genus for halos, giviAx.Haio, 
to that for dark matter, (;max,dm as a function of the 
mean halo mass for each subgroup. The smoothing scale 
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tion, his result reduces to 

ffPT(-.)=-(2^(^) exp(^-^){iJ,(..)- 
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M + S^^'>Hs{i^,) + S^'''>HM'\ a} ,(10) 



where iJ„(f^) are the Hcrmite polynomials and 5^"^ rep- 
resent the skewness parameters, defined as 

3\ _3K!(VV))^ 

4 (7? 



^(2 



9((V;.^-Vi/^)(V^;/^))a3 
4 af 



(11) 



With the relation between i^o- and Vf [equation (5)], the 
genus as a function of is rewritten by 



3PT(l'f) 



1 



(27r)2 ySa 



exp 



' 2 



(5(1) - 5W)i/3(i^f) + (^''^-^^°')ffi (ff )j crj . (12) 

In the case of weak non-Gaussianity induced by the non- 
linear gravitational evolution, the above expression proves 
to be in good agreement with the genus computed from 
the earlier N-body simulation data (Matsubara, Suto 
1996; Matsubara, Yokoyama 1996). We extend those pre- 
vious studies and examine whether or not equations (10) 
and (12) are also applicable to the weak non-Gaussianity 
induced by the halo biasing (in addition to the nonlinear 
gravitational evolution). 

Figure 6 represents a comparison between the simulated 
results for the most strongly biased halo subgroup LL and 
the perturbative model of genus as a function of v^, in 
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Fig. 6. Same as figure 4, but for a comparison of the simulated genus for LL halos (denoted by gualo) with the second-order 
pcrturbative formula [equation (10)] (denoted by Qpt). The skewness parameters in the perturbative formula are calculated in the 
following three ways: (i) a direct computation from the simulated halo density field (solid, labeled by 'sim.'), (ii) fitting of the 
simulated halo genus to the perturbative formula with skewness parameters put as free parameters (dotted, labeled by 'fit'), and 
(iii) direct computation from the <5ijaio field obtained by one-to-one mapping from the dark matter field with the simulated mean 
biasing (dashed, labeled by 'map'). The error-bars, which indicate the sample-to-sample variance [equation(9)], are plotted only for 
the genus with skewness parameters calculated by (i). 



panels (a) and (b) and of i^i in panels (c) and (d). We 
explicitly show the differences between the simulation re- 
sults, gHaio, and the prediction of the perturbative model, 
^PT, divided by the amplitude of the simulated genus, 
ffMAx.Haio, in panel (b) for g^v^) and panel (d) for ^(^'f). 
For the prediction of the perturbative model in equations 
(10) and (12), the r.m.s. fluctuation of the halo density 
field, a = ((5haio^)^^^j and that of its first spatial derivative 
field, (Ti = (|V(5haio^|)^^^, are computed directly from the 
halo density field. We compute the skewness parameters 
S'-°\ S'(^), and S'(^) using three different methods: (i) a 
direct computation from the halo density field (labeled by 
'sim.'), (ii) fitting of our simulated genus to the perturba- 
tive formula as free parameters (labeled by 'fit'), and (iii) a 
direct computation from the Jhaio field (labeled by 'map'), 
based on one-to-one mapping [equation (1)], assuming the 
nonlinear deterministic biasing. For computations of (i) 
and (iii), the gradient and Laplacian values of the den- 
sity fields in equation (11) are evaluated by making the 
corresponding spatial derivative in Fourier space, respec- 
tively. As for (ii) , the range of the fitting is restricted from 
=Max(-l/cr,-2.5) to = 2.5. 

Figure 6 shows that the perturbative formula [equation 
(10)] is valid for the halo genus at Rq > 50/i"^Mpc. The 
approximation by the pcrturbative formula is also good 
in < 1.5 at Rq ~ 30/i^^Mpc, although the large devi- 
ation in > 1.5 might be attributed to the limitation of 
the Edgcworth expansion, which is also seen in the dark 
matter case (Matsubara and Suto 1996). The difficulty 
to distinguish the perturbative prediction using the (5iiaio 



field from other genus implies that the biasing effect on 
the halo genus is well approximated by nonlinear deter- 
ministic biasing. 

Next, we examine whether or not one can reproduce 
the best-fit values for those skewness parameters from the 
independent theoretical models and/or simulation data. 
Figure 7 plots the skewness parameters calculated from 
the above three methods. A weak dependence of the 
skewness parameters on the halo mass is consistent with a 
model prediction by Mo et al. (1997). Fitting method (ii) 
reproduces the skewness parameters well at a smoothing 
scale of 50/i~^Mpc, while the approximation by the per- 
turbative formula breaks down on the scale of 30/i~^Mpc, 
especially for the S^"^^ term. It is almost impossible to re- 
produce skewness parameters accurately at a smoothing 
scale of 100/i^^Mpc, because the statistical error is large 
due to the small value of the genus. 

A good agreement of the skewness parameters from 
the (5haio field again suggests that the approximation by 
the nonlinear deterministic biasing works well in the halo 
genus over the cluster-mass scales. Recalling the fact that 
the nonlinearity in the halo biasing is well described by the 
linear and the quadratic terms shown in section 3, an ac- 
curate prediction of the halo genus is feasible from appro- 
priate modeling of the two biasing coefficients. The reason 
why the stochastic effect of the biasing is nearly negligi- 
ble can be understood from the fact that the stochastic 
components are almost cancelled out by each other in the 
form of skewness parameters (see appendix 1). 
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5. Conclusions 

We evaluate the extent to which halo biasing affects the 
genus for dark halos of cluster-scale mass using the Hubble 
volume simulations. The huge volume of the data enables 
us to perform the most reliable systematic analysis of the 
genus for dark halos and to calculate the sample-to-sample 
error corresponding to the cosmic variance. Through a de- 
tailed comparison of the genus for dark halos with that for 
the mass distribution, we find that the non-Gaussianity 
induced by the biasing is comparable to that by the non- 
linear gravitational evolution. To characterize the biasing 
effect on the genus for dark halos, we extensively apply 
a perturbative formula developed by Matsubara (1994). 
We find that the perturbative formula is indeed very ac- 
curate for the genus of dark halos at a smoothing length 
of i?G ^ 50/i~^Mpc. The skewness parameters obtained 
from a fitting of the perturbative formula consistently re- 
produce the results directly measured from the density 
field. We also find that the nonlinear deterministic bias- 
ing well describes the biasing effect on the genus for dark 
halos. Moreover, the nonlinearity of biasing can be ap- 
proximated by only the linear and quadratic terms. A 
prediction of genus for the dark halo on the cluster-scale 
mass scale is thus possible from an appropriate modeling 
of these two parameters. 

Our results show that the contribution of biasing 
stochasticity to the genus statistics is very small at i? > 
30/i~^Mpc and that the dependence of the halo mass is 
safely negligible within the sampling error-bars. On the 
other hand, both the non-linearity of the biasing and the 
gravitational evolution alter the shape of the genus, de- 
pending on the smoothing scale. In a weakly nonlinear 
regime, one can even discriminate these two effects using 
perturbative theory. In light of this, our formulation, in- 
cluding the nonlinearity of biasing and gravitational evo- 
lution in a previous paper (Hikage et al. 2001), may 
provide a theoretical guide for characterizing the real ob- 
servational data. At the current observational level, how- 
ever, the above biasing effect is very difhcult to detect 
because the biasing effect on the genus is small even in 
a much larger volume size of (1500/i~^Mpc)'^, as shown 
in figure 4. The other observational effects, including the 
redshift distortion and the light-cone effect, are also small 
in the present observational status (Matsubara, Suto 1996; 
Colley et al. 2000). Therefore, the initial Gaussianity can 
be probed directly by the genus from the upcoming galaxy 
cluster catalogs within the sampling errors. 
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Appendix 1. Influence of Biasing Stochasticity in 
a Weakly Nonlinear Regime 

In this appendix, wc discuss the influence of biasing 
stochasticity on the halo genus in a weakly nonlinear 
regime. For this purpose, based on the perturbativc for- 
mula by Matsubara (1994) [equations (10) and (11)], we 
evaluate the skewncss parameters. Suppose that the halo 
density field is split into the deterministic and the residual 
stochastic terms; 



<^halo — <^halo + Shalo, (Al) 

where the mean deterministic term is formally defined by 

mass mass )d5halo, (A2) 

in terms of the conditional PDF of the halo fields at 
a given (5mass- Note that equation (Al) implies that 



(^halo-^(^mass)) 

the relation (Jhai 



for any function F{d, 
oF{S, 



mass 



) because of 



3)) = (^halo-F(^mass))- 

In a weakly nonlinear regime of the dark matter density 
field, the nonlinearity in the biasing may be well approx- 
imated by the following quadratic form: 

b2. 
2 



<5halo('5mass) ~ &l<5mass + (<5n 



(A3) 



with the quantity (Tmm^ being (Jmass )• In this case, the 
variance of the halo and the covariance of the halo and 
dark matter are evaluated as follows: 



2 _ / s- 2 
""hh — \Ohalo 
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\Shalo 



(A4) 

CThm^ = ('5halo'5mass) ~ ^iCTmin^- (A5) 

The above results can be translated to the skewness of the 
halo density field as 
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density field, 5, 
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means the skewness of the mass 
3s'^)/cmm"'- Then, one obtains an 



expression for the normalized skewness, 
density field as 
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where we adopt a measure for the stochasticity of biasing, 
Cscatt (Dekel, Lahav 1999; Taruya, Suto 2000), 



,2 _ 



2 / 2\ 
^mm V'^halo / 



1 (Shalo^ 



(A8) 



In deriving the above results, we did not treat the stochas- 
tic term, Shaio, as a small quantity, which might not be 
guaranteed even in a weakly nonlinear regime. 

Equation (A7) is a natural extension of the frequently 
used expression within the framework of deterministic 
nonlinear biasing (e.g.. Fry, Gaztahaga 1993); 



(0) 



1 



^i?)+3 



bi 



(A9) 



which is only valid in absence of stochasticity, Shaio = 0. 
Compared (A7) with (A9), we readily see the role of 
stochasticity to the skewness parameter as follows: while 
the width of scatter, Cscatt, decreases the skewness am- 
plitude, non-Gaussianity characterized by the third order 
moments, (shaio^'^mass) and (shaio'^) conversely increases 
the overall amplitude. Indeed, this fact illustrates the 
general trend of the skewness parameters. That is, the 
skewness parameters can be schematically divided into 
the contributions from nonlinearity and stochasticity of 
biasing as follows: 



(n) 



bi 



[3rd order stochasticity] 



61 (1 + [2nd order stochasticity]) 



,(A10) 



where the third order stochasticity characterizes the non- 
Gaussianity induced by the stochasticity of biasing and 
the second order stochasticity represents the degree of 
stochasticity. 

Equation (AlO) illustrates the qualitative behavior of 
the normalized skewness for the halo field: (i) it is in- 
versely proportional to the linear bias coefficient 61, (ii) 
the nonlinearity (with 62 > 0) increases the skewness, and 
(iii) the stochasticity decreases the overall amplitude of 
the skewness, but adds a positive term to the purely grav- 
itational term S^. Although it is not clear which of the 
two opposite effects due to the stochasticity is more im- 
portant, the results in figures 6 and 7 show that total 
stochastic effect is not very strong as far as those skewness 
parameters are concerned. This indicates that two oppo- 
site effects due to the biasing stochasticity on the skewness 
parameters [equation (AlO)] are almost cancelled out by 
each other. 



